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Abstract 

We propose a data processing technique that allows searches for a stochastic background of 
gravitational radiation with data from a single detector. Our technique exploits the difference 
between the coherence time of the gravitational wave (GW) signal and that of the instrumental 
noise affecting the measurements. By estimating the auto-correlation function of the data at an off- 
set time that is longer than the coherence time of the noise but shorter than the coherence time of 
the GW signal, we can effectively enhance the power signal-to-noise ratio (SNR) by the square-root 
of the integration time. The resulting SNR is comparable in magnitude to that achievable by cross- 
correlating the data from two co-located and co-aligned detectors whose noises are uncorrelated. 
Our method is general and it can be applied to data from ground- and space-based detectors, as 
well as from pulsar timing experiments. 
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I. INTRODUCTION 



The direct detection of gravitational radiation is one of the most pressing challenges in the 
physics of this century. Predicted by Einstein shortly after formulating his general theory of 
relativity, gravitational waves (GW) will allow us to probe regions of space-time otherwise 
unobservable in the electromagnetic spectrum pQ. Several experimental efforts have been 
underway for awhile, both on the ground and in space [2H9], and only recently kilometer- 
size ground-based interferometers have been able to identify the most stringent upper-limits 
to date for the amplitudes of the radiation expected from several classes of GW sources. 
Although present-generation instruments have not been able yet to unambiguously observe 
GWs, next-generation Earth-based interferometers and pulsar-timing experiments [HI E] , as 
well as newly envisioned space-based detectors [7], are expected to achieve this goal. 

The physical characteristics of the astrophysical sources emitting gravitational radiation 
varies significantly over the observational frequency bands of these detectors. Depending on 
the specific detector design, we will observe GWs emitted during cataclysmic events such 
as supernova explosions and coalescences of binary systems containing neutron stars and 
black-holes, spiraling white-dwarfs binaries in our Galaxy and binary black-holes across the 
Universe, and a stochastic background of astrophysical or cosmological origin pQ. 

The detection of a stochastic background of cosmological origin, in particular, will allow 
us to infer properties about the very early Universe and probe various cosmological models 
predicting the formation mechanism of the background. Since the characteristic strength 
of such stochastic backgrounds is expected to be smaller than the instrumental noise level 
of present and of foreseeable future detectors, it is imperative to identify data processing 
techniques capable of enhancing the likelihood of detection. 

The most general and robust of such techniques relies on cross-correlating the data from 
pairs of detectors operating in coincidence and whose noises are uncorrelated. By virtue of 
observing the same GW signal while being affected by noise sources that are uncorrelated, 
the cross-correlation technique increases the resulting SNR to a stochastic background by 
the square-root power of the integration time. By relying on a network of iV operating 
detectors, the overall sensitivity can be further enhanced by a factor that depends on N 

Although the cross-correlation technique for detecting a background of gravitational ra- 



2 



diation can be applied to ground-based detectors as they are large in number and widely 
separated on Earth (so their noises can be regarded as uncorrelated) , it is presently unthink- 
able to extend it to space-based interferometers as only one might become operational in the 
next decade. Although it might be argued that space-based interferometers such as LISA, 
for example, can generate up to three Time-Delay Interferometric (TDI) combinations [H] 
that could be cross-correlated, it is easy to see that their noises would be correlated making 
the cross-correlation technique ineffective (in simple terms, a TDI combination has some of 
its noises correlated to those affecting the other TDI responses by virtue of sharing with 
them at least "one-arm"). 1 

These considerations imply that a data processing technique for detecting a background 
of gravitational radiation with a single-detector response is required. Although some work in 
this direction has already appeared in the literature for the LISA mission and based on the 
use of the "null-stream" data combination (Midi], further work is still needed. It is within 
this perspective that the present article proposes a new approach to this problem based on 
the use of the auto-correlation function. Our method provides a statistic for enhancing the 
likelihood of detection of a GW stochastic background with a single detector. As we will 
show below, the resulting SNR is comparable to that obtainable by cross-correlating the 
data from two co-located and co-aligned detectors whose noises are assumed uncorrelated. 

The paper is organized as follows. In section [IT] we derive the expression of the SNR 
associated with the auto-correlation function of the data, estimated at an off-set time r that 
is longer than the coherence time of the instrumental noise but shorter than that of the 
stochastic background. After noticing that the SNR can be further enhanced by introducing 
a weighting function, Q(t), in the auto-correlation function, we derive the expression for 



Q that results into an optimal SNR. Finally in section |III| we present our comments and 
conclusions, and emphasize that our auto-correlation technique can be applied not only to 
the data from a forthcoming space-based interferometers, but also to those from a network 
of ground-based interferometers and pulsar timing. 



1 Although it is possible to construct, within the TDI space, three TDI combinations whose noises are un- 
correlated [T2] , it is straightforward to show that also their responses to an isotropic stochastic background 
of gravitational radiation would be uncorrelated. 
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II. THE AUTO-CORRELATION FUNCTION 



Let us denote with R(t) the time-series of the detector output data, which we will assume 
to include a random process, n(t), associated with the noise and also a random process 
associated with a GW stochastic background, h(t), in the following form 

R(t) = h(t) + n(t) . (1) 

In what follows we will assume the GW stochastic background to be an unpolarized, station- 
ary, and Gaussian random process with zero mean. As a consequence of these assumptions 
such a background is characterized by a one-sided power spectral density, P h (\f\) defined by 
the following expression [TD1 US] 

(h(f)h*(f')) = l$(f-f')Ph(\f\), (2) 

where the symbol ~ represents the operation of Fourier transform, the * symbol denotes the 
usual operation of complex conjugation, and the angle-brackets, (), denote the operation of 
ensemble average of the random process. We will further assume n to be stationary, Gaussian 
distributed with zero- mean, and characterized by a one-sided power spectral density, P n (\f\), 
defined as follows 

(n{f) n*{f)) = l -5{f - f) P n {\f\) . (3) 

Let us denote with r gw and r n the coherence time of the GW stochastic background and 
of the instrumental noise respectively, and let us also assume r gw > r n . This assumption is 
rather general since, by being true for white noise and any colored GW background, it must 
remain true in general as the data can always be pre- whitened 2 . 

The sample auto-correlation function, S(t), of the data, R(t), can be written in the 
following form 

f T/2 

S{t) = / R(t) R(t +r)dt, (4) 

J-T/2 

where T is the time interval in seconds over which the auto-correlation is computed. The 
interval of interest for our method is when r n < r < r gw ; in this region the instrumental 

2 Although the case when both the noise and the stochastic background spectra are white over the entire 
observational band cannot be tackled by our technique, theoretical models for likely backgrounds of cos- 
mological and astrophysical origins characterize them with "colored" spectra that are in general different 
from those associated with the noises of the detector [T] 
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FIG. 1: Schematic illustration of the method. The upper panel shows the power spectrum of 
(assumed- white) instrumental noise plus a low- frequency excess ("red") signal power. The lower 
panel shows the auto-correlation function of the spectrum in the upper panel, expanded around 
zero correlation. The white component decorrelates rapidly from unity, while the red process gives 
rise to long-range correlations. The detection of these correlations competes with estimation error 
statistics of the correlation function, with an SNR increasing like T 1//2 , as described in the text. 

noise has decorrelated, leaving the correlation function of the signal. Figure [T] shows the 
situation schematically. This signal correlation then competes with the estimation error 
statistics of the correlation function for detection, as discussed below. Very similarly to the 
cross-correlation technique discussed in [TOj HE] we can rely on S(t) for building a decision 
rule for the detection of a GW stochastic background. 

From the statistical properties of both the GW stochastic background and the noise, and 
by taking into account the central-limit theorem (for r gw << T, S(t) is equal to the sum of 
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many uncorrelated products) it follows that S(r) is also a Gaussian random process [15]. In 
order to fully characterize it we can estimate its mean and variance, which must be functions 
of only the spectra of the two random processes h(t) and n(t). 

By taking the ensemble average of S and S* 2 over many noise and GW signal realizations 
and for T » r gw , after some long but straightforward calculations we get the following 
expressions for the mean and variance of S(t) (see Appendix [A] for details) 3 

HS = (S) = - J J h {\f\)e 2mfT df , (5) 

rp POO 

4 = (S 2 ) - (S) 2 = ^ J cos 2 (27r/r) + P n (|/|)] 2 df . (6) 

Eqs.(|5j [6]) can then be used for deriving the following expression for the Signal-to- Noise 
Ratio (SNR auto ) 



a s V 2 



SNR auto = ^ = \l- J '°° — ■ (7) 



/_ 00 00 cos 2 (2 7 r/r) [P h (\f\) + P n (|/|)] 2 df 



In order to get a quantitative understanding of the above expression for the SNR of 
the auto-correlation at lag r, we may compare it against the SNR associated with the 
cross-correlation (computed also at time-lag r) of the data from two co-aligned and co- 
located detectors whose noises are assumed to be uncorrelated (see [HI] for details). It is 
straightforward to obtain the following expression for the cross-correlation signal-to-noise 
ratio SNR 2 d 

T n^l/De 2 -^/ 

Too K(l/D cos 2 (2vr/r) + \P h (\f\) (Pm(|/|) + + I^CI/D^d/l)] 1 ^ 

(8) 

where we have denoted with Pj n (|/|) the one-sided power spectral density of the noise in 
detector i = 1,2. In order to compare the two expressions given in Eqs. (jTj, [s]) , let us assume 
the two noise power spectral densities P in to be equal to each other and to a white spectrum 
having spectral level equal to A per unit bandwidth. By taking the ratio between SNR 2 d 



3 It should be noticed that Eq. (JgJ) can also be derived from the equation for the estimation error of the 
auto-correlation function given in Jenkins and Watts |17j (eq. 5.3.20), under the approximation that 
T » r gw , and that the width of the square of the auto-correlation function is small compared with T. 
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and SNR auto we get 



SNR 



2d 



JZcos^nfr) [P h (\f\) + A] 2 df 



1/2 



SNR, 



auto 



IZo [PhdfD cos*(2tt/t) + P h (\f\) A + \A*\ 



1/2 



(9) 



Since current estimates for the likely amplitudes of cosmological gravitational wave back- 
grounds indicate them to be significantly smaller than the instrumental noise levels of 
current and future detectors, if we take P^ « Am. Eq. ^ in this limit we find 
SNR 2 d/ SNR auto -> l. 4 



A. The optimal filter 

Although the physical reasons underlying the selection of the time interval within which 
the time-lag r should be selected are quite clear, they do not allow us to uniquely identify a 
value for it. In particular, since the number of possible time lags we could use for estimating 
the auto-correlation function can be in general rather large, it is natural to look for the 
optimal way of combining all possible auto-correlation functions. In analogy with the cross- 
correlation technique [TU] in which an optimal filter Q(t,t') maximizing the SNR could be 
identified, here we show that also for the auto-correlation technique there exists such a 
filter. Since the stochastic background and the instrumental noise have been assumed to 
be stationary, we introduce a time-invariant filter function and write the following linear 
combination of all possible sample auto-correlation functions lags 

pT/2 p+oo 

S= / R(t)R{t + r)Q(r) dt dr , (10) 

J -T/2 J-oo 

where the auto-correlation is computed over the time interval T, and the filter function 
Q(t) at the time-lag r is non-zero in the interval r n < \t\ < r gw and zero everywhere else (a 
property that we have already taken into account in Eq. (10) by extending the r-integration 
over the entire real axis). 

From the statistical properties of both the GW stochastic background and the noise, and 
by virtue of the central-limit theorem, it follows that S is also a Gaussian random process 
that can be fully characterized by estimating its mean fig and variance, <r|. Although 

It is easy to show that, also in the limit of "high signal-to-noise ratio", Ph >> A, S N R,2d / S N R auto —> 1. 
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their derivations are quite long, they are straightforward and result into the two following 
expressions 

rp p+QO 

^=2] p h(\f\) Q(f) d f > ( n ) 

rp i>+oo 

4 = 2 J Wl/I) + p «(l/l)] 2 IW)I 2 # • ( 12 ) 



By simple inspection of Eq. ( 12 ) , and following [10J , it is convenient to define the following 



operation of inner product between two arbitrary complex functions, say A, and B 

/+oo 
A(f)B*(f)[Ph(\f\) + Pn(\f\)] 2 df . (13) 
-oo 

With this newly defined inner product, n$ an d a s can be rewritten in the following form 

w= l^iafki i) ; ^=! (e '«' (14) 

while the squared signal-to-noise ratio of S, SNR 2 is equal to 

S ^-4"2 (Q,Q) ■ (15) 



From a simple geometrical interpretation of the inner product defined through Eq. (13) 



from Eq. (15) it is easy to see that the maximum of the SNR of S is achieved by choosing 



the filter function Q to be equal to 

Q(/) = A mmfkm ' (16) 

where A can be any arbitrary real number. The above expression for the filter Q implies the 
following maximum SNR achievable by the auto-correlation technique 

SNRlax = 2 [p.d/D+^d/i)] 2 df ■ (17) 

III. DISCUSSIONS AND CONCLUSIONS 



The main result of our work has been to show that the auto-correlation function of the 
data measured by a single detector of gravitational radiation can be used for searching for a 
stochastic gravitational wave background. This is possible when there is a difference between 
the spectral properties of the GW background and of the instrumental noise. By estimating 
the auto-correlation function of the data at a time-lag at which the auto-correlation of the 
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noise has decayed to a value smaller than that of the auto-correlation function of the GW 
background, one can enhance the resulting power signal-to-noise ratio by the square-root of 
the integration time. 

We have assessed the effectiveness of our technique by comparing the expression of its 
SNR against that achievable by cross-correlating the data from two co-located and co-aligned 
detectors whose noises are assumed to be uncorrelated. We found the two expressions to 
coincide in the expected limit of the noise power spectrum larger than that of the GW 
background. In order to further enhance the SNR associated with the auto-correlation 
technique we have derived the expression of an optimal filter function that reflects any prior 
knowledge we might have about the spectrum of GW background. 

The auto-correlation technique can be used when searching for a stochastic background 
of gravitational radiation with (i) networks of ground-based interferometer as it will further 
enhance the likelihood of detection achievable with the cross-correlation method, (ii) with a 
single space-based interferometer, and (iii) with pulsar timing data. 
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Appendix A: Derivation of the mean and variance of S(t) 

In what follows we derive the expressions of the mean and variance of the auto-correlation 
function S(t) (Eqs. ([5j [6]) given in Section ([ii])). We will assume the instrumental noise to 
be characterized by a random process n(t) that is Gaussian-distributed with zero- mean, 
and has one-sided power-spectral density P n (\f\). We will also take the GW stochastic 
background to be a Gaussian random process of zero-mean, unpolarized, and stationary. 
As a consequence of these assumptions such a background is uniquely characterized by a 
one-sided power spectral density, Ph(\f\) as given by Eq. §2§. 

From the definitions of the data, R(t), and of the auto-correlation function, S(t), calcu- 
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lated at a time-lag r such that r n < t < r gw , we have 

•T/2 



S(t)= [ [h(t) + n(t)] [h{t + r) + n(t + t)] dt , (Al) 

J -T/2 

where T is the time interval in seconds over which the auto-correlation function is computed 
and it is assumed to be longer than r gw . After taking the ensemble average of both sides of 



Eq. (Al ), we get 



r-T/2 

(S(t)) = I [(h(t)h(t + r)> + (h(t)n(t + r)> + (n(t)h(t + r)> + (n(t)n(t + r))] dt . 

-T/2 

(A2) 



Note that the second and third terms on the right-hand-side of Eq. (A2) are equal to 
zero because the random process associated with the stochastic GW background and that 
associated with the instrumental noise are uncorrelated. The fourth term instead can also 
be regarded to be zero because is the auto-correlation function of the noise estimated at an 
off-set time longer than its coherence time. From these considerations we find the following 
expression for the mean of S 

rT/2 rp p+OO 

^ = (S(r)> = / [{h(t)h(t + r)>] dt =- P h (\f\)e 2mfT df . (A3) 

J -T/2 Z J-oo 

In order to derive the expression for the variance of S (Eq. ^) we need to calculate (S 2 ). 
Since S 2 (t) can be written as 

r-T/2 r T/2 

i2, 



S 2 (r) = / / [h(t)h(t + t) + h(t)n(t + t) + n(t)h(t + t) + n(t)n(t + t)] 

J -T/2 J -T/2 

x [h(t')h{t' + t) + h(t')n(t' + r) + n(t')h(t' + t) + n(t')n(t' + r)] dt dt' , (A4) 



by taking the ensemble average of both sides of Eq. (A4) only nine of the sixteen terms 
present are different from zero. This follows from the Gaussian nature of both random 
processes and because terms such as (x{t)y{t')z(t + r)w(t' + r)) (where x,y,z,w can be 
either h or n) can be expressed in the following form [15] 



(x(t)y(t')z(t + r)w(t' + r)> = (x(t)y(t'))(z(t + r)w{t' + r)> + (x(t)z(t + r)){y{t')w{t' + r)> 

+ (x(t)w(t' + T))(y(t')z(t + r)) (A5) 
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From the above considerations we get the following expression for (S 2 (t)) 

r T/2 rT/2 

(S 2 (t)) = / / [(h(t)h(t + r)){h(t')h(t' + r)> + (h(t)h(t'))(h(t + T)h(t' + t)) 

J-T/2 J-T/2 

+ (h(t)h(f + r))(h(t')h(t + r)} + (h{t)h{t')){n{t + r)n(t' + r)> 
+ (/i(t)/i(t' + r)){n{t + r)n(t')} + + r))(n(t)n(t' + r)) 

+ (h(t + T)h(t' + r)) (n(t)n(0) + (n(t)n(t')){n(t + r)n(t' + r)> 
+ (n(t)n(t / + r))(n(t / )^(t + r))] dt dt' , (A6) 

which can further be rewritten in the following form 

rp2 r+oo r+oo rp r+oo 

(S 2 (r)) = — / P h (\f\)P h (\f\) e 2 ^ f ~ f ' )T df df' + - P»(\f\) df 

^ J—oc J — oo ^ «/ — oo 

m /»+oo on /*+QC 

+ Z / ^ 2 (l/l)e 4m/r df + - \ P h (\f\)Pn(\f\) df 

^ J —oo ^ J —oo 



-oo ^ </ — oo 

-OO rp r+oo 



rp r+OO rp r+OO 

+ - / P h (\f\)P n (\f\)e^ df+- P h (\f\)P n (\f\)e-^ df 

^ J — oo ^ «/ — OO 

rp /»+oo m /*+co 

+ j Ph{\f\)Pn{\f\)df + - P 2 n {\f\)df 
^ J —00 ^ J —00 

rp r + OO 

+ T / ^(l/l)e 4m/T d/ • (A7) 
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00 



If we now square both sides of Eq. (A3) and subtract it from Eq. (A7) to obtain the 
expression for the variance cr|, after some rearrangement of the remaining terms we get the 
following result 

rp r+oo 

4=2/ cos 2 (27r/r)[P,(|/|) + P n (|/|)] 2 df . (A8) 
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